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Abstract 


y output vector 


The application of the eigensystem synthesis 
technique to place the closed-loop eigenvalues 
and shape the closed-loop eigenvectors has not 
been practical for active flutter suppression, 
primarily because of the availability of only 
one control surface (aileron) for flutter sup- 
pression. The oblique-wing aircraft, because of 
its configuration, provides two independent sur- 
faces (left and right ailerons), making the appli- 
cation of eigensystem synthesis practical. This 
paper presents the application of eigensystem 
synthesis using output feedback for the design 
of an active flutter suppression system for an 
oblique-wing aircraft. The results obtained are 
compared with those obtained by linear quadratic 
Gaussian techniques. 

Symbol s 

A pi ant matrix 

B control matrix 


output matrix 
Identity matrix 
eigensystem gain matrix 


partitioned matrix of specified 
components of (A-ji - A)"*B 


L' complex conjugate transpose of L 

M partitioned matrix of unspecified 

components of (Ail - A) _1 B 

u input vector 

V eigenvector matrix 


* 

Ai , Ai closed-loop eigenvalue and its 

conjugate 

a_ minimum singular value of the return 

difference matrix 

Introduc tion 

The U.S. Navy and NASA are currently involved 
in the design and development of an unsymmetric- 
skew-wing aircraft capable of 65° wing sweep and 
flight at Mach 1.6. Such a unique configuration 
exhibits aeroelastic behavior distinctly dif- 
ferent than that of straight, swept-back, or 
swept-forward wings and has a potential for poor 
modal response characteristics. The most serious 
result of such characteristics can be flutter, an 
unstable motion caused by an interaction between 
structural vibrations and aerodynamic forces. 1 
Active suppression of aerodynamic wing flutter 
can result in substantial weight savings and 
increases in performance compared with passive 
methods such as increased structural stiffness 
and mass balancing. 

The synthesis of modal control systems for 
unsymmetric aircraft requires considerably more 
states than are necessary for symmetric config- 
urations because all degrees of freedom must be 
adequately represented in a single formulation. 
The model representing the aircraft must include 
rigid-body modes, flexible modes, unsteady aero- 
dynamics, actuators, and gust states. Control 
laws have been formulated for active flutter con- 
trol using the standard linear quadratic Gaussian 
(LQG) procedure; 2-4 the synthesized optimal feed- 
back control laws are of the same order as the 
aircraft plant. Practical izing the control law 
requires order reduction by techniques such as 
transfer function matching, modal truncation, and 
residual ization.^ - ? 


Vi closed-loop eigenvector 

v?, v? achievable and desired eigenvectors, 

respectively 

* 

Vi complex conjugate of Vi 

x state vector 

*Associate Professor, Electrical Engineering 
Department. Member AIAA. 

**Aerospace Engineer. 


To investigate modal control synthesis strate- 
gies for an oblique-wing configuration, a generic 
skew-wing aircraft model was developed for 45° 
wing skew at a flight condition of Mach 0.70 at 
10,000 ft altitude. At this flight condition the 
aircraft has an unstable flutter mode. An active 
impl ementabl e flutter control law was developed^ 
using the LQG design technique and modal resid- 
ualization to reduce the order of the controller. 
However, this method increased the order of 
response of the closed-loop system as compared 
with that of the open-loop system and degraded 
modal control characteristics. 



Eigensystem synthesis procedures are suitable 
for flight control system design because they do 
not increase the order of the system. Also, the 
difficulty in incorporating specifications such 
as damping, frequency, and decoupling within a 
quadratic performance index makes the eigensystem 
synthesis procedure a promising design alterna- 
tive. The performance specifications can be 
interpreted in terms of desired closed-loop eigen- 
values and eigenvectors. Moore^ and others have 
shown how feedback can be used to place closed- 
loop eigenvalues and shape closed-loop eigenvec- 
tors. References 10 to 12 successfully demon- 
strate the use of an eigenstructure assignment 
procedure for aircraft control system design. 

The eigensystem synthesis technique using 
output feedback to place closed-loop eigenvalues 
and shape closed-loop eigenvectors has not been 
used for active flutter suppression, primarily 
because of the availability of only one control 
surface for flutter suppression, making eigen- 
vector shaping impractical. The oblique-wing 
aircraft, because of its configuration, provides 
two independent control surfaces (left and right 
ailerons) and makes the application of eigensystem 
synthesis practical. 

This paper presents the application of 
eigensystem synthesis to the design of an 
active flutter suppression system for a generic 
model of an oblique-wing aircraft. The results 
obtained are compared with those obtained by 
LQG techniques. 

Flutter Suppression Usi ng E igenstructur e 

The generic oblique-wing aircraft aerostruc- 
tural model used in the system synthesis process 
is a simple beam representation of the fuselage 
and wing. The aerodynamic model was developed by 
superimposing aero panels over the beams, as shown 
in Fig. 1. The aircraft modal characteristics 
were developed using NASTRAN analysis. At the 
sweep configuration (45°) and flight condition 
(Mach 0.70, 10,000 ft) selected, the unaugmented 
aircraft has a flutter mode characterized as 
primarily wing bending. Because the intent of 
this paper is to demonstrate a design synthesis 
process, the model order was reduced considerably; 
the final model contained a rigid-body (primarily 
pitch) mode along with three elastic modes. The 
model reduction process did not significantly 
affect the flutter mode characteristics. Details 
of the aircraft model formulation are given in 
Ref. 8. 

The aircraft design model includes the 
linearized form of the unsteady aerodynamics, 
actuators, and gust model dynamics and can be 
represented as 


x = Ax + Bu 

(1) 

<< 

1! 

O 

X 

(2) 


where x is a 24 * 1 state vector , y is a 5 * 1 out- 
put vector, u is a 2 x l input vector, and A, B, 
and C are the plant, control, and output matrices, 
respectively, of suitable dimensions. The 
24 states include the rigid-body mode, flexible 
mode deflections, flexible mode rates, unsteady 


aerodynamic states, actuator deflection and rate 
states, and wind gust states. Eight states result 
from the retained structural modes, eight from 
the two-lag-term set of approximated unsteady 
aerodynamics, six from the two actuators, and 
two from the gust model. The five outputs con- 
sist of pitch angle, pitch rate, and three accel- 
erations: center of gravity and right and left 

wingtips. The two inputs are right and left 
aileron deflections. 

For the system under consideration, the 
following conditions hold;H 

1. A maximum of five closed-loop eigenvalues 
can be assigned arbitrarily with the stipulation 
that if A-j is a complex closed-loop eigenvalue, 

★ 

its complex conjugate A-j must also be a closed- 
loop eigenvalue. 

2. A maximum of five eigenvectors can be 
altered. If the shape of a complex eigenvector 

* 

vi is altered, its complex conjugate vi must 
be altered in the same way. 

3. For each eigenvector whose shape is 
altered, a maximum of two eigenvector elements 
can be chosen arbitrarily. 

4. Achievable eigenvectors must lie in the 

subspace spanned by the columns of (Ail - A) _1 8 of 
dimension two, which is the number of independent 

control variables. A desired eigenvector v*? will 

in general not reside in the prescribed subspace 
and cannot be achieved. The optimal achievable 

eigenvector v(j is obtained by the orthogonal 

projection of v^ onto the subspace spanned by 

(AjI - a) _ 1 B. 

Given the system of equations (1) and (2) and 
assuming an output feedback, the control input 
u is given by 

u = Ky (3) 

where K is the gain matrix of dimension 2 x 5 . 

For the closed-loop system, the following rela- 
tionship holds: 

(A - BKC)vi = A-j vi (4) 

where vi is a closed-loop eigenvector and Aj a 
closed-loop eigenvalue; or 

(Ail - A)vi = -BKCvi = BMi (5) 

where 

Mi = -KCvi (6) 

In general , the desired eigenvector v^ does not 
reside in the achievable subspace, and an approxi- 
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mate solution by methods of orthogonal projection 
can be obtained. 13 

If 

Li = (Ail - A) _1 B (7) 

then the achievable eigenvector v? is given by 

A - LitLiLiT'LiV? (8) 

and the gain K is given by 

K = -M(CV) t (9) 

where the superscript t denotes pseudoinverse, L^ 
is the complex conjugate tranpose of Li, and 


M = [Mi 

M 2 ... 

M n 3 

(10) 

V = [vj 

a 

v 2 ... 

a n 

v n ] 

(ID 


Result s 

The control law synthesized by the method 
outlined was applied to the design of an active 
flutter suppression controller for an oblique-wing 
aircraft. A generic 45° wing-skew structural 
model was developed to simulate flutter at a sub- 
sonic flight condition of Mach 0.70 and an alti- 
tude of 10,000 ft . 

Table 1 gives the eigenvalues of the open-loop 
aircraft model. The unstable eigenvalue pair at 
this flight condition (0.5 ± 14 . 37i ) represents 
primarily wing bending. 

The design objective was to stabilize the 
aircraft without exceeding the specified root- 
mean-square (rms) control activity so that sat- 
uration would not occur. Based on actuator 
limitations, the rms deflection of the aileron 
was limited to 5° and the deflection rate to 
30 deg/sec. In addition to stabilizing the air- 
craft with low surface activity, it is required 
that the controller be robust. The controller 
considered here is multi-input, multi-output: the 

right and left wing control surfaces are indepen- 
dent of each other because of the unsymmetric 
nature of the aircraft. Robustness of the multi- 
loop control system is evaluated by using the 
singular values of the return difference matrix. 34 

The process of selecting desired closed-loop 
eigenvalues and eigenvectors for the case under 
consideration presents a problem. The selection 
of desired eigenvalues and eigenvectors normally 
is based on engineering judgment and requires a 
clear insight into physical aspects of the plant 
being controlled. Most of the physical insight 
for the flutter problem is lost in the process of 
model reduction from the original large dimensions 
to the 24th-order model developed for the control 
system synthesis process. 


A time-consuming and arbitrary method for the 
selection of eigenvectors did lead to stabiliza- 
tion of the aircraft. Another method, using 
closed-loop eigenvalues and eigenvectors obtained 
through LQG design assuming full -state feedback, 
was used to obtain a stable closed-loop system. 

The closed-loop eigenvalues and eigenvectors 
obtained by this method were selected as the 
desired eigenvalues and eigenvectors for the 
eigensystem synthesis process. 

Table 2 gives the desired locations of the 
five closed-loop eigenvalues as well as the 
desired eigenvectors. Five eigenvalues were 
located, including the unstable flutter mode, at 
the desired location. The achieved eigenvectors 
are given in Table 2. The feedback gain K was 
evaluated based on the achieved eigenvalues and 
eigenvectors. The values of K and the closed- 
loop eigenvalues of the aircraft are also indi- 
cated in Table 2. 

Table 3 gives the comparative description of 
the rms values of response to gust input for both 
the LQG design8 (seventh-order controller) and the 
controller designed using eigensystem synthesis. 
The eigensystem synthesis rms values for the sur- 
face deflections and rates are within the pre- 
scribed range and compare favorably with those 
obtained by the other method. 

The controller developed by this method (K is 
a 2 x 5 matrix) is extremely simple to implement. 
In comparison, the LQG design technique uses a 
full -order controller (the order being the same as 
that of the plant). Even a reduced-order 
controller is difficult to implement and increases 
the order of the system. 

However, this eigensystem synthesis approach 
compromises robustness. Stability robustness of 
multi-input, multi-output feedback control systems 
is characterized by the minimum singular value of 
the return difference matrix of the plant input or 
output. 14 Figure 2 shows the plots of the minimum 
singular values a_ for the controller developed in 
this paper and the full-state feedback controller 
developed using the LQG technique. 8 A degradation 
in robustness is evident from the plot. However, 
the reduced-order controller is not as robust as 
the one designed using the eigensystem procedure, 
as is evident from Fig. 2. 

Concl usions 

An impl ementabl e flutter controller for a 
45°-skew oblique-wing aircraft mathematical model 
was designed using the eigensystem synthesis tech- 
nique. The controller does not increase the order 
of the system and is extremely simple to imple- 
ment, whereas the LQG design technique uses a 
full-order controller, is difficult to implement, 
and increases the order of the system. Work is in 
progress to improve stability margins by using 
constrained optimization techniques to shape the 
singular value spectrum. 35 A standard performance 
index is minimized while trying to satisfy minimum 
singular value constraints at the plant input or 
output , or both . 
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Table 1 Open-loop 
eigenvalues 


0.0000 


O.OOOOi 

-0.4187 

+ 

O.OOOOi 

-0.4229 

+ 

O.OOOOi 

-4.2075 

+ 

O.OOOOi 

-0.1612 

+ 

5.0036i 

-0.1612 

- 

5 . 0036i 

-1.9397 

+ 

14. 0551 i 

-1.9397 

- 

14.0551i 

0.5011 

+ 

14. 3649i 

0.5011 

- 

14.36491 

-20.0000 

+ 

O.OOOOi 

-20.0000 


O.OOOOi 

-30.8526 

+ 

O.OOOOi 

-35.1822 

+ 

0.5223i 

-35.1822 

- 

0. 5223i 

-37.1170 

+ 

3. 564 7 i 

-37.1170 

- 

3. 564 7 i 

-39.6927 

+ 

0.6080i 

-39.6927 

- 

0.6080i 

-41.3346 

+ 

O.OOOOi 

-36.4000 

+ 

37. 1354i 

-36.4000 

- 

37.13 54i 

-36.4000 

+ 

37 . 1354i 

-36.4000 

- 

37.1354i 
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Table 2 Eigensystem variables 


Desired eigenvectors (from LQG method) Desired eigenvalues 


3.9999 

0.0438 

+ 

0.0049i 

0.0438 

- 

0.0049i 

0.0015 

+ 

0.0007i 

0.0015 

- 

0.0007i 

-0.0003 

+ 

O.OOOOi 

0.0000 

-0.0317 

+ 

O.OOlSi 

-0.0317 

- 

0.0018i 

-0.0245 

+ 

0.00891 

-0.0245 

- 

0.0089i 

-0.2355 

- 

5.0020i 

0.0000 

0.0011 

+ 

0.00011 

0.0011 

- 

O.OOOli 

0.0023 

+ 

0.0004i 

0.0023 

- 

0.0004i 

-0.2355 

+ 

5.0020i 

0.0000 

0.0047 

+ 

0.00011 

0.0047 

- 

O.OOOli 

-0.0004 

+ 

0.0004i 

-0.0004 

- 

0.0004i 

-0.5041 

- 

14. 3657 i 

-0.0010 

-0.0349 

+ 

0 . 21 7 7i 

-0.0349 

- 

0. 2177i 

-0.0104 

+ 

0.0216i 

-0.0104 

- 

0.02161 

-0.5041 

+ 

14. 3657i 

0.0000 

-0.0017 

- 

0.1 588i 

-0.0017 

+ 

0.1588i 

-0.1154 

- 

0.3565i 

-0.1154 

+ 

0.3565i 




0.0000 

-0.0009 

+ 

0.0057i 

-0.0009 

- 

0.0057i 

-0.0075 

+ 

0.0330i 

-0.0075 

- 

0.0330i 




0.0000 

-0.0014 

+ 

0.0236i 

-0.0014 

- 

0.0236i 

-0.0054 

- 

0.00651 

-0.0054 

+ 

0.0065i 




0.0746 

2.8445 

+ 

16.26161 

2.8445 

- 

16.26161 

-6.5976 

- 

0.962H 

-6.5976 

+ 

0.9621i 




0.0612 

1.9139 

+ 

5.3855i 

1.9139 

- 

5.3855i 

2.4618 

+ 

6.5556i 

2.4618 

- 

6 . 5 5 56 i 




-0.0247 

-7.3179 

- 

74.96821 

-7.3179 

+ 

74.96821 

2.4855 

- 

3.2472i 

2.4855 

+ 

3.2472i 




0.0829 

0.9809 

- 

29.58321 

0.9809 

+ 

29. 5832i 

10.8184 

- 

32.7421i 

10.8184 

+ 

32.742H 




-0.0699 

-2.6257 

- 

15.7088i 

-2.6257 

+ 

15. 7088i 

7.1162 

+ 

3.8695i 

7.1162 

- 

3 . 8695 i 




-0.0578 

-1.6892 

_ 

3.96591 

-1.6892 

+ 

3.9659i 

-2.0765 

- 

5 . 652 5i 

-2.0765 

+ 

5.6525i 




0.0015 

6.1769 

+ 

73.08021 

6.1769 

- 

73.0802i 

-2.6786 

- 

0.0123i 

-2.6786 

+ 

0.01231 




-0.1024 

-1.9195 

+ 

26.9120i 

-1.9195 

- 

26. 9120i 

■12.8369 

+ 

31.10181 

-12.8369 

- 

31.10181 




0.0000 

0.0000 

+ 

O.OOOli 

0.0000 

- 

O.OOOli 

0.0004 

+ 

0.0003i 

0.0004 

- 

0.0003i 




0 . 0000 

-0.0003 

+ 

O.OOOli 

-0.0003 

- 

O.OOOli 

-0.0050 

+ 

0.0055i 

-0.0050 

- 

0.0055i 




0.0000 

-0.0006 

- 

0 . 001 6i 

-0.0006 

+ 

0.0016i 

-0.0766 

- 

0. 0740i 

-0.0766 

+ 

0.0740i 




0.0000 

0.0000 

+ 

O.OOOOi 

0.0000 

- 

O.OOOOi 

0.0009 

- 

O.OOOli 

0.0009 

+ 

O.OOOli 




0.0000 

0.0000 

+ 

O.OOOOi 

0.0000 

- 

O.OOOOi 

0.0016 

+ 

0.0130i 

0.0016 

- 

0.0130i 




0.0000 

-0.0001 

- 

O.OOOli 

-0.0001 

+ 

O.OOOli 

-0.1878 

+ 

0.01 59i 

-0.1878 

- 

0.01 59i 




0.0000 

0.0000 

+ 

O.OOOOi 

0.0000 

- 

O.OOOOi 

0.0000 

- 

O.OOOOi 

0.0000 

+ 

O.OOOOi 




0.0000 

0.0000 

+ 

O.OOOOi 

0.0000 

- 

O.OOOOi 

0.0000 


O.OOOOi 

0.0000 

+ 

O.OOOOi 








Achieved eigenvectors 






Achieved closed-loop 
eigenval ues 

3.9999 

0.0438 

+ 

0.0049i 

0.0438 

_ 

0.0049i 

0.0015 

+ 

0.0007i 

0.0015 

_ 

0.0007i 

-0.0003 

+ 

O.OOOOi 3 

0.0000 

-0.0317 

+ 

0 . 001 8i 

-0.0317 

- 

0.0018i 

-0.0245 

+ 

0.0089i 

-0.0245 

- 

0.0089i 

-0.4187 

+ 

O.OOOOi 

0.0000 

0.0011 

+ 

O.OOOli 

0.0011 

- 

O.OOOli 

0.0023 

+ 

0.0004i 

0.0023 

- 

0.0004i 

-0.4229 

+ 

O.OOOOi 

0.0000 

0.0047 

+ 

O.OOOli 

0.0047 

- 

O.OOOli 

-0.0004 

+ 

0.0004i 

-0.0004 

- 

0.0004i 

-3.6775 

+ 

O.OOOOi 

-0.0010 

-0.0349 

+ 

0 . 21 7 7i 

-0.0349 

- 

0 . 21 77i 

-0.0104 

+ 

0.0216i 

-0.0104 

- 

0 . 021 6i 

-0.2355 

+ 

5. 0020i 3 

0.0000 

-0.0017 

- 

0.1 588i 

-0.0017 

+ 

0. 1588i 

-0.1154 

- 

0.35651 

-0.1154 

+ 

0. 3565i 

-0.2355 

- 

5.00201 3 

0.0000 

-0.0009 

+ 

0.0057i 

-0.0009 

- 

0.0057i 

-0.0075 

+ 

0.0330i 

-0.0075 

- 

0.0330i 

-1.3818 

+ 

12. 944 5i 

0.0000 

-0.0014 

+ 

0.0236i 

-0.0014 

- 

0.0236i 

-0.0054 

- 

0.0065i 

-0.0054 

+ 

0.0065i 

-1.3818 

- 

12. 944 5i 

0.0746 

2.8445 

+ 

16. 261 6i 

2.8445 

- 

16.26161 

-6.5976 

- 

0 . 962 li 

-6.5976 

+ 

0.96211 

-0.5041 

+ 

14. 36571 a 

0.0612 

1.9139 

+ 

5.3855i 

1.9139 

- 

5.3855i 

2.4618 

+ 

6.5556i 

2.4618 

- 

6.5556i 

-0.5041 

- 

14. 3657 i a 

-0.0247 

-7.3179 

- 

74.9682i 

-7.3179 

+ 

74. 9682i 

2.4855 

- 

3 . 2 4 7 2 i 

2.4855 

+ 

3 . 2472 i 

-23.1705 

+ 

O.OOOOi 

0.0829 

0.9809 

- 

29. 5832i 

0.9809 

+ 

29.5832i 

10.8184 

- 

32 . 7421i 

10.8184 

+ 

32. 7421i 

-21.2819 

+ 

13. 991 3i 

-0.0699 

-2.6257 

- 

15. 7088i 

-2.6257 

+ 

15. 7088i 

7.1162 

+ 

3.8695i 

7.1162 

- 

3.8695i 

-21.2819 

- 

13. 991 3i 

-0.0578 

-1.6892 

- 

3.9659i 

-1.6892 

+ 

3.9659i 

-2.0765 

- 
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-2.0765 

+ 

5.6525i 

-34.5807 

+ 

1 . 041 Oi 

0.0015 

6.1769 

+ 
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6.1769 

- 
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- 
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+ 
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-34.5807 

- 
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-1.9195 
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- 
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+ 
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- 
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+ 
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+ 
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- 
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+ 
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0.0004 
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+ 
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0.0000 

- 

O.OOOOi 
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,0000 0.0004 

-0.0006 

-0.0001 

0.0002 
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,0000 0.0009 

-0.0010 

-0.0001 

-0.0001 







asee desired eigenvalues. 
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Table 3 Root-mean-square responses at flutter conditions 



Right 

wi ng 

Left 

wi ng 


6 , deg 6 

, deg/sec 

6 , deg 

6 , deg/sec 

Full -state feedback 

1.58 

7.91 

0.28 

3.58 

Full-order controller 
with Kalman estimator 

1.58 

12.58 

0.41 

5.64 

Full-order controller 
with robust Kalman 
estimator 

2.06 

11.21 

0.41 

5.02 

Reduced-order 
control 1 er 

1.51 

10.04 

0.41 

4.90 

Eigensystem 

0.62 

4.81 

0.69 

8.46 


synthesi s 



Fig. 1 Generic model (aero panels and node points). Vertical fin 
shown in X-X plane. 



Fig. 2 Minimum singular values of the return dif- 
ference matrix. 
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